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We present results of a simulation of 2 flavour QCD on a 16 3 x 4 lattice using p4- 
improved staggered fermions with bare quark mass m/T = 0.4. Derivatives of the ther- 
modynamic grand canonical partition function Z{V,T, fi u , fi d ) with respect to chemical 
potentials n u> d for different quark flavours are calculated up to sixth order, enabling es- 
timates of the pressure and the quark number density as well as the chiral condensate 
and various susceptibilities as functions of fi Ut d via Taylor series expansion. Results are 
compared to high temperature perturbation theory as well as a hadron resonance gas 
model. We also analyze baryon as well as isospin fluctuations and discuss the relation to 
the chiral critical point in the QCD phase diagram. We moreover discuss the dependence 
of the heavy quark free energy on the chemical potential. 



1. Introduction 



It is important to study QCD at high temperature and high density by numerical 
simulations of lattice QCD. In particular, studies of the equation of state (EoS) can 
provide basic input for the analysis of the experimental signatures for QGP formation, 
e.g. the EoS will control the properties of any hydrodynamic expansion. We extend 
studies of the EoS to non-zero baryon number density. 

Simulation at non-zero baryon density is known to be difficult; however recent studies 
have found that a Taylor expansion with respect to chemical potential [i u ^ is an efficient 
technique to investigate the low density regime, interesting for heavy- ion collisions [QJI2]- 
In the calculation of the Taylor expansion coefficients, i.e. calculation of the derivatives 
at Hu,d = 0, the technical difficulty for non-zero \i u ^ does not arise and a quantitative 
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Figure 1. The ratios of Taylor expansion coefficients for pressure (quark number suscep- 
tibility), isospin susceptibility and chiral condensate. 



study becomes possible within the error by truncation of higher order terms. Since ther- 
modynamic quantities can be defined by derivatives of the partition function, e.g. chiral 
condensate (4>ip) = (T/V)(d\n Z/dm), quark number density n q = (T/V)(d\n Z/dfi q ), 
quark number susceptibility Xq — (T/V)(d 2 \iaZ/d{i 2 ), and isospin susceptibility xi — 
(T /V)(d 2 In Z/dfij), where \i q = (fi u + (J>d) /2 and fj,j = — /id)/2, the calculation of first 
derivatives yields basic QCD thermodynamics observables, and the calculation of higher 
derivatives is a natural extension. 

In this study, we evaluate the Taylor expansion coefficients of thermodynamic quantities 
and compare these results to the predictions from the perturbation theory in the high 
temperature limit and from the hadron resonance gas model in the low temperature 
hadron phase. Also fluctuations of quark number, isospin number and electric charge are 
discussed near the transition temperature T c . They are estimated by Xqi Xi an d charge 
susceptibility xc = X<?/36 + Xi/4, and are related to the critical endpoint expected at 
non-zero Moreover we study the Taylor expansion coefficients of chiral condensate 
and heavy quark free energy. The details of simulations are given in [Ej. 



2. Quark gluon gas and hadron resonance gas 

We define the expansion coefficients by p/T A = (\nZ)/(VT 3 )_ = J2™=o C n{.^q/T) n , 
Xq[Xi]/T 2 = En= 2 n(n - l)c n [4](^/T)- 2 and (^>/T 3 = J^dl+fa/T)" for N = 0. 
We expect the equation of state to approach that of a free quark-gluon gas (Stefan- 
Boltzmann (SB) gas) in the high temperature limit. The coefficients in the SB limit 
for 2 flavour QCD with fij = are well known as c 2 = c\ = 1,C4 = c{ = l/(27r 2 ) and 
c n = = for n > 6. Moreover the sign of Cq is negative starting at 0(g 3 ) in pertur- 
bation theory. On the other hand, in the low temperature phase QCD is well described 
by a hadron resonance gas model in which the pressure is obtained by summing over 
the contributions from all resonance states of hadrons. In this model, the contribution 
to p/T 4 from each individual hadron which has baryon number Bi is in proportion to 
exp(3Bifi q /T), hence the pressure can be written as p/T 4 = G(T) + F{T) cosh(3/i g /T), 
where G(T) and F(T) are the mesonic and baryonic components of p/T 4 at fi q = 0. 
Similarly, we obtain Xq /T 2 = 9F(T) cosh(3/i g /T), X i/T 2 = G 7 (T) + F 7 (T) cosh(3/^/T) 
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Figure 2. The quark number susceptibility Xq/T 2 (left) and isospin susceptibility Xi/T 2 
(right) for various \i q jT . T is T c at \i q = 0. 

and (i/j^/T 3 = G^{T) + F^(T) cosh(3/x 9 /T). Here, the mesonic component for Xq is 
zero because mesons have baryon number zero. Therefore, C4/C2 = 3/4, cq/c^ = 3/10, 
cl/ c { = 3/10 and cf/cp = 3/4 in the region where the hadron resonance gas provides a 
good approximation. 

We investigate these coefficients. The results for c n+ 2/c n , c^ +2 /c^ and c^j_ 2 / c n are 
shown in Fig. ^ We find that these results are consistent with the prediction from the 
hadron resonance gas model for T/T c < 0.96 and approach the SB values, i.e. c 4 /c 2 = 
c l/ c 2 = l/^ 71 " 2 , c 6 /c 4 = Cq/cI = 0, in the high temperature limit. Also c 6 is negative 
at high temperature as expected in the perturbation theory. These results suggest that 
the models of free quark-gluon gas and hadron gas seem to explain the behaviour of 
thermodynamic quantities well except in the narrow regime near T c . 

3. Susceptibilities at non-zero \i q 

Next, we calculate quark number and isospin susceptibilities in a range of < fi q /T < 1. 
The data connected by solid lines in Fig. El are obtained by Xq/T 2 = 2c2 + 12c,i(fi q /T) 2 + 
30c6(/i g /T) 4 and the corresponding equation for xi- Dashed lines are the results from 
0(fi 2 ) expansion. Since the statistical error of cq is still large near T c , the location and 
height of the peak are less accurately determined, and also the error from the truncation 
of higher order terms of u. q /T seems to be visible for large n q /T. However, as seen in 
Fig. [H cq changes its sign at T c . This means the peak position of Xq moves left, which is 
corresponding to the change of T c as a function of \i q . T c (u, q /T = 1)/T c (fi q /T = 0) in [ 
HH is about 0.93. Moreover, Xq increases with higher orders of the expansion for T < T c 
which confirms the existence of a peak. This suggests the presence of a critical endpoint in 
the (T,n q ) plane. On the other hand, xi i n Fig. El does not show any singular behaviour. 
This is consistent with the sigma model prediction that only isosinglet degrees of freedom 
become massless at the critical endpoint. 
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Figure 3. (a) The ratios ml/ml and mf/ml as functions of T/T c and (b) the singlet 
screening mass at non-zero \i q . 



4. Screening mass at non-zero ji q 

Finally we want to discuss the free energy of a static quark anti-quark pair at finite 
temperature and density. We extract the singlet free energy F^q and colour averaged free 
energy Fqq by the Polyakov loop correlation functions. For T > T c , the free energy is ex- 
pected to be exponentially screened at large distances, AF^ 1 (r, T, fi q ) = F^q (oo, T, fi q ) — 
FgQ(r,T,fj, q ) ~ e -«i av ' 1 ( 7 >9) r, - \Ye calculate the Taylor expansion coefficients of the De- 
bye screening mass m 3 *' 1 in terms of \i q jT at \i q = 0, m 3 *' 1 = J2^=o m n V,1 (/ x g/^) n - We 
plot the data of m\/m\ and jm\ in Fig. Efa). We find m\ = (l/2)m| v . This is ex- 
pected from perturbation theory, which suggests that the leading order contribution to 
the colour singlet free energy is given by one gluon exchange while the colour averaged free 
energy is dominated by two gluon exchange. Moreover the perturbative Debye screen- 
ing mass mo is given by mD(T,fi q ) = rriDfi(T)^Jl + 3Nf/[(2N c + Nf)TT 2 ](^ q /T) 2 , with 

mD,o{T) = g(T)T^J (2N C + iV/j/6. Here iV c and Nf are the number of colour and flavour 
respectively. The solid line in Fig. Efa) is the perturbative prediction for m\/mQ. The 
ratio m\/m\ is found to be consistent with perturbation theory for T > 2T C . In Fig. E^b) 
we show the \i q - dependence of the singlet screening mass m 1 (fi q ,T)/T for a small values 
of Hq/T, where only contributions from and m\ are included. Further details of this 
study are given in 
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